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Aims. We use a sample of 83 core-dominated active galactic nuclei (AGN) selected from the MOJAVE (Monitoring of Jets in 
AGN with VLBA Experiments) radio-flux-limited sample and detected with the Fermi Large Area Telescope (LAT) to study 
the relations between non-simultaneous radio, optical, and 7-ray measurements. 
■ Methods. We perform a multi-band statistical analysis to investigate the relations between the emissions in different bands and 
I reproduce these relations by modeling of the spectral energy distributions of blazars. 
.* . Results. There is a significant correlation between the 7-ray luminosity and the optical nuclear and radio (15 GHz) luminosities 
of blazars. We report a well defined positive correlation between the 7-ray luminosity and the radio-optical loudness for 
^—^ J , quasars and BL Lacertae type objects (BL Lacs). A strong positive correlation is found between the radio luminosity and the 
I ■ 7-ray-optical loudness for quasars, while a negative correlation between the optical luminosity and the 7-ray-radio loudness is 
• • present for BL Lacs. Modeling of these correlations with a simple leptonic jet model for blazars indicates that variations of the 
accretion disk luminosity (and hence the jet power) is able to reproduce the trends observed in most of the correlations. To 
reproduce all observed correlations, variations of several parameters, such as the accretion power, jet viewing angle, Lorentz 
factor, and magnetic field of the jet, are required. 
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1. Introduction 

Blazars are observed over a wide range of the electromag- 
netic spectrum from radio to very-high-energy (VHE) 7- 
rays. Synchrotron emission, which is believed to be domi- 
nant in the radio band, originates in the relativistic bipo- 
lar outflows which on milliarcsecond scales (parsec scales) 
appear as moving supcrluminal radio knots, which can be 



* Table [4] is only available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5). 



probed by very-long-baseline inte rferometry (VLBI) im- 
ages (e.g.. iKellermann et al. 1998 ). 

One qu arter of a ll AGN of the radio-flux-limited 



MOJAVE (jLister et al.l l2009al ) sample objects were de- 
tected by the Fermi LAT during the first three months 
of operation ([Lister et al. 2009b). For this sample of 31 
blazars a significant correlation between radio (VLBA0) 
flux density at 15 GHz and 7-ray photon flux was re- 
ported for quasi-simultaneous radio - 7-ray measurements 
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( Kovalev et al. 20091) . They suggested that the radio and 
7-ray emissions are produced in the cores of parsec-scale 
jets. A positive correlation between radio f lux density at 
8.4 G Hz and 7-ray flux was reported by iLinford et al. 
( 201lh for 50 flat-spectrum quasars selected from the 



VLB A Imaging and Polarimetry Survey (VIPS). The cor- 
relation appeared to be stronger for radio bright quasars 
and the differences between the 7-ray loud and quiet flat- 
spectrum quasars can be explained by Doppler boost- 
ing. Strong correlation between the 7-ray flux (above 100 
MeV) and the 20 GHz flux density was found for 134 flat- 
spectrum quasars and BL Lacs which were identified from 
cross correlation of the 1FGL with the 20 GHz Australia 
Telescope Compact Array radio s urvey catalogue (Mahon i 



et al. 2010, Ghirlanda et al. 2011). lPushkarev et all (|2010l) 
studied 187 bl azars selected from the First Fermi LAT 
catalog (1FGL; lAbdo etah 2010b ) and found that on av- 
erage the 7-ray flares precede the radio flares at 15 GHz 
with intrinsic time delays of 1.2 months. This result sug- 
gests that the origin of the bulk of the the 7-ray photons 
is located within the 15 GHz VLBI core. High-frequency 
VLBI (43 GHz) and multi- waveband monitoring of blazars 
indicate that a single superluminal radio knot can cause 
a number of 7-ray flares, and that very rapid variability 
of 7-ray emission co mes near the location of t he 43 GHz 
cores of blazars jets ( Marscher fc Jorstad 2010f ). Recently, 



Leon-Tavares et al.l (|2011fi have shown that the strongest 
7-ray flares follow the mm flares in about 70 days. They 
interpret this as strong 7-ray events being produced at a 
distance of several parsecs downstream of the mm-wave 
core. 

The MOJAVE blazars show a tight connection be- 
tween radio VLBA and optical emission. Long-term radio 
- optical monitoring of radio galaxies showed that a vari- 
able optical emission on scales from a few months to a few 
years is generated in the innermost 0.4 pc region of the ra - 
dio jet ( Arshakian et al. 2010a ; Leon-Tavares et al. 2010l ). 
A significant positive correlation at milliarcseconds scales 
was reported between the optical nuclear emission and (i) 
the radio emission of the unresolved cores at 15 GHz in 
superluminal quas ars, and (ii) the extended jet emission in 
BL Lac objects by I Arshakian et aL I (l2010bh . Those results 
suggested that the radio and optical emissions are beamed 
and originate in the innermost part of the sub-parsec-scale 
jet in quasars. 

Extreme optical variability was a clear signature 
for some blazars detected at 7-ray energies during the 
CGRO/EGRET era (Fichtel et al. 1994). Taking advan- 
tage of the improved gamma sensitivity of Fermi /LAT, 
recent multiwavelength variability studies on individ- 
ual sources have confirmed the tight connection among 
the 7-ray and var iable optical emission (e.g., 3C273 in 



Abdo et a l. 2010a). 



So far there was no attempt to compare optical contin- 
uum properties among 7-ray blazars. In this manuscript 
we study multi-band correlations to investigate the single 
production mechanism for continuum emission in differ- 
ent wavebands. Taking advantage of the unprecedented 



sensitivity provided by the 1FGL, we analyze relations 
between optical nuclear emission, radio emission on scales 
of milliarcseconds, and 7-ray emission of more than 80 
MOJAVE AGN that had high confidence Fermi /LAT 
source associatio ns after its fi rst 11 months of operation. 
As discussed in Abdo et ail ()2010bl ). it is possible that 
some blazars may lack a Fermi /LAT association due their 
location in a particularly confused location (numerous 7- 
ray sources, or high diffuse background near the galac- 
tic plane). Since the MOJAVE sample includes only very 
radio-bright sources and avoids the galactic plane region, 
we will refer for convenience in this paper to the non-LAT 
associated MOJAVE sources as 'non-detected'. An inter- 
pretation of the correlations and the main physical process 
driving these correlations is discussed in the framework of 
a steady-state jet model. 

The sample of blazars is defined in Sect. 2. In Sect. 3, 
the Kendall's r statistical test is used to examine the cor- 
relations between emission characteristics defined in radio, 
optical, and 7-ray regimes. We interpret the implications 
of the jet model in Sect. 4 and summarize our results in 
Sect. 5. 

Throughout the paper a flat cosmology with H = 71 
km s _1 Mpc -1 , O m = 0.27, and £1a = 0.73 is adopted. 

2. The sample of AGN detected with Fermi/LAT 

We study a statistically complete sample of 135 core- 
dominated AGN compiled from the MOJA VE sample and 



the VLBA 15 GHz monitoring surve y (Kcllcr mann et al 



19981 : IZensus et ail l2002t iLister et all l2009al) . The com- 
plete sample, MOJAVE-1, is limited by the total VLBA 
(mas-scale) flux density at 15 GHz: it is greater than 1.5 
Jy for sources in the Northern hemisphere (S > 0) and > 2 
Jy in the Southern hemisphere (—20 < S < 0). The sam- 
ple comprises quasars, BL Lacs, radio galaxies and some 
unidentified sources having typically a core-jet structure 
on milliarcsecond scales. At high redshifts, the sample has 
a selection bias towards BL Lacs and quasars with fast and 
bright one-sided jets oriented closer to the line of sight, 
while at low redshifts (z < 0.1) there is a selection bias 
towards slow and intrinsically bright jets of radio galaxies 
(for more details see Lister & Homan 2005, Arshakian et 
al. 2010b). 

The 1FGL catalogue includes 83 sources from the 
MOJAVE-1 sample (hereafter, M1-1FGL sample): 59 
quasars, 19 BL Lacs, 3 radio galaxies, and 2 unidentified 
sources. Most of the blazars not detected at 7-ray ener- 
gies come from the weaker radio portion of the MOJAVE-1 
sample (Fig. [IJ indicating that the 7-ray selected sample 
is biased towards radio strong blazars. The latter also are 
strong 7-ray blazars since there is a strong positive corre- 
lation between radio and 7-ray fluxes of blazars (Kovalev 
et al. 2009). 

For these 83 7-ray selected sources, optical nuclear 
fluxes a nd redshift measureme nts were available for 76 
objects ( Arshakian et alJfeoiObf) : 58 quasars, 15 BL Lacs, 
and 3 radio galaxies (see Tabled]). A procedure to correct 
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Fig. 1 Distribution of radio VLB A flux densities at 15 GHz 
for 135 MOJAVE-1 AGN (full line) and AGN detected 
with Fermi /L AT (filled blue histogram) for all blazars, 
quasars, BL Lacs, radio galaxies, and optically unidenti- 
fied sources. 



th e photometric fluxes for a stellar emission is discussed 
. Arshakian et al. I (|2010bh . Typical errors associated with 



optical nuclear fluxes are between 15% and 25%. 

We derive the optical nuclear luminosities (L opt ), 
total VLBA luminosities (£vlba), and the rest frame 
radio-optical loudness (the ratio of the fC-corrected ra- 
dio and optical fluxes, R — <S' VLBA /S' opt = <SVlba(1 + 
z)~ ai+a ° / S op t, where a r — is the spectral index of the 
total VLBA component, a a is the optical spectral index 
adopted to be —0.5). We use the optica l nuclear fluxes at 



5100 A given in I Arshakian et al.l (|2010bf ) and the total ra- 



dio flux densities at 15 GHz measured (during the period 
from 1994 to 2003 as part of the VLBA 15 GHz monitor- 
ing survey) at the epochs at which the unresolved compo- 
nent flux density has its maximum value. 7-ray luminosi- 
ties w ere computed using the Eq. (1) in IGhisellini et al 
(|2009h . 



And 



, 2 S-yjvx,^) 
(1 + z) 1 ~ a 'r 



(1) 
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Fig. 2 Distributions of optical nuclear fluxes (top panel) 
and luminosities (bottom panel) at 5100 A for 122 AGN 
from the MOJAVE-1 sample (full line). Among these, 76 
AGN are detected as the gamma-ray source (filled his- 
togram) with Fermi /L AT during the 1FGL period. 



where S^iyx, V2) is the energy flux between v\ =0.1 GeV 
and v-2 — 100 GeV from the 1FGL catalogue, a 7 is the 7- 
ray energy index between the frequencies v\ and i>2, and 
d,L is the luminosity distance. Note that the 7-ray, opti- 
cal, and radio luminosities of the MOJAVE-1 blazars are 
estimated from non-simultaneous observations and their 
typical uncertainties are between ~ (10 — 15)% in 7-ray 
band , ~ (15 — 25)% in optical, and ^ 5% (jKovalev et al 



20051 1 in radio band. The rest frame 7-r ay-radio loudness 



75- 



VLBA 



is defined as G r — S 1 /Sy LBA — S 7 (l 
and the 7-ray-optical loudness is G = 
z)- a -< +a °/S opt . 

We use the non-parametric Kendall's r test to analyze 
correlations betw een independent variable s and the partial 
Kendall's r test (jAkritas fc Siebertl Il996h to account for 
the common dependence on redshift in the correlations 
between luminosities. Throughout the paper we assume a 
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Fig. 3 Optical nuclear luminosity against 7-ray luminos- 
ity (top panel), and radio VLBA (15 GHz) luminosity 
against 7-ray luminosity (bottom panel) for AGN from 
the M1-1FGL sample. Labels in the top panel denote the 
corresponding population. 

correlation to be significant if the chance probability for a 
null hypothesis is P < 0.05. 

3. Relations between 7-ray, optical, and radio 
properties of superluminal AGN 

Detection rate of j-ray blazars. The relative number of 
blazars detected in 7-rays increases with increasing radio 
flux, and most quasars not detected in 7-rays are radio 
weak (Fig. [IJ . This is explained by a tight correlation be- 
tween radio and 7-ray emission of quasars (Kovalev et al. 
2009). Most of BL Lacs and radio bright galaxies are de- 
tected in 7-ray band. 

The distributions of optical nuclear fluxes for detected 
and not detected 7-ray blazars have similar ranges (Fig. [2] 
top panel). The Kolmogorov-Smirnov (K-S) statistical test 
shows that the null hypothesis that these distributions 
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Fig. 4 7-ray luminosity (L 7 ) versus radio-optical loudness 
(R = S VLBA /S opt ) for the M1-1FGL quasars and BL Lacs. 
The dotted and dashed lines represent the ordinary least- 
square ht to the data for quasars and BL Lacs respectively. 



are drawn from the same parent population can not be 
rejected (P — 0.15; confidence level of 85 %). We conclude 
that among radio-selected blazars, the detection rate of 7- 
ray AGN does not depend on their optical nuclear fluxes, 
and most of undetected 7-ray blazars are in fact quite 
optically luminous (Fig. [2J bott om panel) 
Optical 



Radio emission. lArshakian et al.l (|2010bh 
found a positive correlation in the L pt ~ -^vlba plane for 
AGN from the MOJAVE-1 sample. They concluded that 
the correlation is due to the population of quasars and that 
the optical emission is non-thermal and generated in the 
parsec-scale jet. This is supported by st udies of individ 



ual ra d io galaxies, 3C 390.3 and 3 C 120 (jArshakian et al 



2010al iLeon-Tavares et al.l l20ld) for which a link be 



tween optical continuum variability and kinematics of the 
parsec-scale jet was found. It was interpreted in terms of 
optical continuum flares generated at subparsec-scales in 
the innermost part of a relativistic jet rather than in the 
accretion disk. 

Gamma-ray - Optical emission. No correlation is 
found for BL Lacs, while a significant positive correlation 
between L 1 and L opt is found for quasars (99%; Fig. [31 
Table [IJ suggesting a single production mechanism for 7- 

ray and optical nuclear emission. 

(|2009h 



Kovalev et al 



Gamma-ray - Radio emission. 
reported a correlation between 7-ray and radio VLBA 
(8 GHz) emission for a sample of ~ 30 AGN. This corre- 
lation holds at high confidence level (> 99.9%; Fig. [3]) for 
non-simultaneous measurements at 15 GHz and a larger 
sample (M1-1FGL) of AGN. This suggests that the pow- 
ers averaged over the long time scales are correlated and, 
hence, the Doppler-factors of the parsec-scale jet in the 
gamma and radio domains are n ot changing subst a ntially 
on a timescale of a few years. IPushkarev et al.l (|2010h 
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Table 1 Kendall's r correlation analysis between emission characteristics of AGN from the M1-1FGL sample. Al and 
A2 are the independent variables for which the Kendalls r correlation analysis is performed, and A3 is the control 
variable (if exists then the partial Kendalls r correlation analysis is applied to Al and A2), r is the correlation 
coefficient, and P is the probability of a chance correlation. The correlations are considered to be significant if the 
chance probability P < 0.05 (or confidence level > 95%). 
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found that the 7-ray emission leads the radio emission of 
the parsec-scale jet at 15 GHz with a time delay of a few 
months. They interpreted the observed time lag as a re- 
sult of synchrotron opacity in the jet: the radio and 7-ray 
emissions can be ge nerated in the same re gion by per- 
turbations of the jet ( Marscher Sz Gear 19851 ) and become 
observable with some time delay due to the opacity effects. 
If this scenario is correct, the variable optical emission is 
also generated in the perturbation moving upstream in the 
jet, and one should expect that the optical emission leads 
the radio emission and is also delayed with respect to the 
7-ray emission. 

The Lvlba ~ L 1 correlation is visibly tighter (0.39 dex 
scatter) than the L opt - L 1 one (0.54 dex scatter; Fig. [3]). 
The remarkably tight relation between Lvlba and I/ 7 , 
logivLBA = (0.83 ± 0.03) log L 7 + (5.96 ± 0.26), over a 
range of about five orders of magnitudes, can be used to 
estimate the 7-ray flux of not detected 7-ray blazars from 
their VLBA flux density. One may think t hat the corrected 



optica l nuclear emission of some AGN ([Arshakian et al 



2010bJl might be contaminated by a contribution of a stel- 
lar component, thus causing the large dispersion in the 
L op t ~ L-y diagram. This contamination should be stronger 
in radio galaxies and weaker in quasars and BL Lacs for 
which the contribution of optical nuclear emission is dom- 
inant. The large dispersion in the L opt - L 1 plane can be 
due to stronger variability in the optical regime than in 
the radio, and/or a wider range of Doppler factors in the 
optical regime compared to the range of Doppler factors 
of the jet at the location where the 15 GHz emission is 
produced, if the bulk of the optical emission is generated 
in the relativistic jet and it is Doppler boosted. 

We report a significant positive correlation (> 99 %) 
between L 7 and radio-optical loudness R for quasars and 
BL Lacs (see Fig. [Hand Table[T]). For the combined sample 



the L~. 



R' 



0.95±0.11 



BL Lacs 



R 



0.79±0.27 



for quasars L 7 oc /j - 44 ± 018 ) an d f or 
The slope derived for quasars 
is shallower than for BL Lacs. To test the robustness of 
this result, we excluded one outlier quasar (logi? = 1.31 
and logL 7 = 45.63) and repeated the fitting of data. The 
slope for quasars increases, 0.58 ±0.17, and coincides with 



the slope of BL Lac objects within the error range. We 
conclude that the slopes for quasars and BL Lacs are 
not statistically different. A noticeable difference in radio- 
optical loudness is found between populations of quasars 
and BL Lacs: logi? = 0.12 for quasars and logi? = —0.81 
for BL Lacs. 

Gamma-ray loudness. We define the rest frame 7- 
optical loudness (G ) and 7-radio loudness (G r ) as G = 
Sry/S^t and G r = S^/ S vhBA , respectively. We find that 
the radio luminosity and 7-optical loudness are signifi- 
cantly correlated only for quasars (see Table [T] and Fig. [5j 
left panel), £ V lba oc G°- 45±012 . 

The correlation in the L opt - G r plane (Fig. [5J right 
panel) is marginally significant only for BL Lacs (Tabled]). 
A large amount of scatter about the fitted line and small 
sampling of BL Lacs account for the large error of the 
power index (L opt oc G-° S1±0A8 ). 

Logarithm of the median 7-optical loudness is signif- 
icantly different for quasars and BL Lacs (2.1 and 1.2), 
while there is no difference between their medians of 
7-radio loudness (2.1 and 2). There are no high-energy 
peaked BL Lacs in the MO J AVE- 1 sample. The latter 
have significantly different SED properties than the 
radio-selected MOJAVE-1 BL Lacs (Lister et al. 2011). 
The differences found between quasars and BL Lacs are 
also supported in recent studies (e.g., Ghisellini et al. 
2009, Sambruna et al. 2010, Tornikoski et al. 2010, Gupta 
et al. 2011) suggesting the presence of different physical 
conditions along the jet in quasars and BL Lacs. This and 
other insights into the relationship between the 7-ray, 
radio, and optical emission will be pursued in more detail 
in further studies. 

Relations found between L 7 and R, Lvlba and G , 
and L op t and G r are due to variations of the SED of 
blazars, which can be driven by changes of the Doppler 
factor, magnetic field, and intrinsic power of the jets. In 
the next section, we attempt to identify the main physical 
parameters of the jet responsible for the observed correla- 
tions. 
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Fig. 5 Le/t panel. Radio VLBA luminosity against 7-optical loudness (G a = S„/S opt ) for the M1-1FGL quasars and 
BL Lacs. Right panel. Optical luminosity against 7-radio loudness (G r = ^/Svlba)- The dotted and dashed lines 
represent the ordinary least-square fit to the data for quasars and BL Lacs respectively. 



4. Interpretation and model implications 

The broadband emission of blazars is generally understood 
as originating from ultrarelativistic particles in the rela- 
tivistic jet, which is oriented at a small angle 6> b s ~ 1/r 
(where F is the bulk Lorentz factor of the material moving 
along the jet) with respect to our line of sight. The low- 
frequency continuum emission from radio to optical - UV 
(or even X-rays, in high-frequency-peaked BL Lac objects) 
frequencies is generally accepted to be synchrotron emis- 
sion from relativistic electrons in the jet. The high-energy 
(X-ray-7-ray) emission can be either inverse Compton 
emission from the same electrons (leptonic models), or 
emission resulting from the production of secondaries in 
photo-pair and photo-pion production reactions of ultra- 
relativistic hadrons in the jet (hadroni c models) . For a 
recent review of blazar models, see, e.g., iBottcherl (j2010l ). 



We will discuss a possible interpretation of the cor- 
relations described above, in the framework of a steady- 
state, one-zone leptonic model. Specifically, we adopt the 
semi-anal ytical description of a gen eric leptonic blazar 
model of iBottcher fc Dermerl ( 20021 ). In this model, a 
power-law distribution of relativistic electrons, Q(j) = 
Qo 7~ 9 i?(7; 71, 72) with low- and high-energy cut-offs 71 
and 72 is persistently injected into a spherical emission re- 
gion of radius Rb, which moves with constant speed /?r c, 
corresponding to bulk Lorentz factor T, along the jet. In 
order to reduce the free parameter space, we choose an 
observing angle of 6* b s = 1/r so that the Doppler fac- 
tor D = (T [1 — Pr cos#obs])~ = T. A tangled magnetic 
field B is present in the emission region. Modeling of the 
SEDs of many BL Lac objects indicates that a model in- 
cluding only synchrotron and synchrotron-self-Compton 
(SSC), is often sufficient t o describe those SEDs (e.g 



2000l : Krawczvnski et al. 2002 ). Therefore, we adopt a 
pure SSC model in this study to represent the SEDs of 
high-frequency-peaked blazars (BL Lac objects), which 
is fully determined th rough the model paramete rs listed 
above. In the model o f lBottcher fc Dermerl (j2002l) . the in- 
jected electron distribution is followed in an iterative pro- 
cess to equilibrium between ongoing particle injection and 
radiative cooling and escape. Fig. [6] shows a characteristic 
sequence of BL Lac SEDs resulting from a variation of the 
jet power Lj , while keeping the magnetic field at equipar- 
tition with the equilibrium electron distribution. Typical 
parameters for one of those SED calculations are listed in 
Table 03 

For low-frequency-peaked BL Lacs (LBLs) and flat- 
spectrum radio quasars (FSRQs), it is generally accepted 
that leptonic models require external photon sources 
for Compton scattering (EC = Extern al Compton 



1997; Sambruna et al. 


Madeiski et al. 


1999 


Hart man et al. 


2001). 



emis si on (e.g. 



Dermer et al 



Il997t iMukheriee et all 
Bottcher fc Blooinl 



1999; 



2000; 



Mastichiadis fc Kirkl fl997l IPian et al.l Il998t IPetrv et al 



In the case of quasars, the 
most prominent external radiation field might be the 
quasi-isotropic radiation field of the broad-line region 
(BLR), as long as the emission region is located within 
the BLR. We model the BLR as a spherical shell at a 
distance i?BLR from the central engine, reprocessing a 
fraction tblr of the accretion disk luminosity, Lrj. In this 
case, the jet power is parameterized as a fraction fj of 
the accretion disk luminosity, i.e., Lj — fj Lrj. As in the 
case of BL Lac objects, we assume the persistent injection 
of a power-law distribution of ultrarelativistic electrons 
into the emission region, and let the distribution evolve 
to an equilibrium between injection, radiative cooling, 
and escape. A typical set of parameters producing an 
FSRQ-like SED is listed in Table El A sequence of FSRQ 
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Fig. 6 Sequence of model SEDs for BL Lac objects, vary- 
ing the jet power, while adjusting the magnetic field to 
equipartition with the relativistic electrons in the jet. The 
curves are labelled by the value of Lj. The black dotted 
and dashed lines indicate the individual radiation compo- 
nents adding up to the the total spectrum with the highest 
Lj = 10 43 erg s _1 . This sequence yields the points shown 
as blue open squares in Figs. 151 -[TU1 . 

SEDs, resulting from a change in the accretion disk 
power (and hence also the jet power), while adjusting the 
magnetic field to be in equipartition with the equilibrium 
electron distribution, is shown in Fig. [7J 

We tested whether variations of individual parameters 
for FSRQs and BL Lacs were able to reproduce the ob- 
served statistically significant (P < 0.05) correlations be- 
tween R and L 7 , G and Lvlba for quasars, and between 
G r and L opt for BL Lacs. In particular, we tested whether 
these correlations can arise from a change in the following 
fundamental parameters: (a) Accretion disk (and hence, 
jet) power Ljj (Lj), (b) Accretion disk (and jet) power, 
while keeping the magnetic field in equipartition with the 
equilibrium electron distribution (es = us/u e — 1), (c) 
Accretion disk (and jet) power with cb — 1, adjusting 

1/2 

the BLR radius as -Rblr cx L^ , (d) Lorentz factor T 
(and, hence, the Doppler factor D), (e) viewing angle 6* b s 
(also changing the Doppler factor D), and (f) magnetic 
field B. The resulting theoretical correlations for FSRQs 
are shown by the filled symbols in Figs. |8j - [101 while 
open symbols show representative sequences for BL Lac 
objects. In all cases, we ran one sequence per base param- 
eter. Figs. 151-lTTIlshow representative results. The straight 
lines in these figures indicate the best power-law fits to 
the observed correlations. 

When comparing the theoretical correlations with the 
observed ones, one needs to keep in mind that our 
model only takes into account radio emission from the 
high-energy emission region in the sub-pc scale core. 
Realistically, we expect a substantial contribution to the 
observed radio emission to arise on larger (> pc) scales. 
Therefore, our model radio fluxes should be considered as 




[Hz] 



Fig. 7 Sequence of model SEDs for quasars, varying the 
accretion disk luminosity Ljj (and hence, the jet power 
Lj = fjLjj), while adjusting the magnetic field to equipar- 
tition with the relativistic electrons in the jet. The curves 
are labelled by the value of Lp ■ The black dotted, dashed, 
dot-dashed, and dot-dot-dashed lines indicate the individ- 
ual radiation components adding up to the total spectrum 
with the highest Lo = 10 46 erg s _1 . This sequence yields 
the points shown as blue filled squares in Figs. 151- 1TU1 



Table 2 Typical leptonic jet model parameters for a 
BL Lac object. 



Quantity 


Symbol 


Value 


Minimum injection electron Lorentz factor 


71 


10 4 


Maximum injection electron Lorentz factor 


72 


10 6 


Electron injection spectral index 


q 


2.5 


Magnetic field 


D 


2 G 


Blob radius 


Rb 


3 X 10 15 cm 


Bulk Lorentz factor 


r 


15 


Observing angle 


Sobs 


3.82° 


Jet power 




10 42 erg s -1 



Table 3 Typical leptonic jet model parameters for a 
quasar. 



Quantity 


Symbol 


Value 


Minimum injection electron Lorentz factor 


71 


300 


Maximum injection electron Lorentz factor 


72 


10 5 


Electron injection spectral index 


q 


3.1 


Magnetic field 


D 


0.661 G 


Blob radius 




10 17 cm 


Bulk Lorentz factor 


r 


15 


Observing angle 




3.82° 


Accretion disk luminosity 




10 45 erg s" 


Jet power fraction 


fj 


0.05 


Reprocessing Optical Depth of BLR 


T BLR 


0.1 


Effective radius of BLR. 




10 1S cm 



1 



lower limits. In particular, additional radio emission from 
larger scales would shift all model curves to the right in 
Fig. [8] and to the left in Fig. [9] Considering this potential 
offset, we identified a set of parameters with which the ob- 
served R vs. Lj (Fig. [5J, G r vs. L op t (Fig. O, and G vs. 
£vlba (Fig. [TU|) correlations can be reasonably well repro- 
duced through a variation of Lp (and Lj), while keeping 
the magnetic field in equipartition with the electron dis- 
tribution. Those sequences are marked by the blue dashed 
line (square symbols) in Figs. |5J - QUI and the sequences 
of SEDs are illustrated in Figs. [5] and [7] for typical pa- 
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Fig. 8 Model correlations of radio-optical loudness R vs. gamma-ray luminosity L 1 along sequences varying individual 
parameters. Each sequence represents the effect of a variation of one individual base parameter. The legend in the top 
left indicates the parameters varied along the sequences. Filled symbols indicate SEDs characteristic of FSRQs, while 
open symbols indicate sequences of SEDs characteristic of BL Lacs. The long-dashed and dot-dot-dashed lines mark 
the power-law correlations that best fit the observational data. The sequences marked in blue dashed line (squares), 
showing reasonable agreement with the best-fit correlations, correspond to the sequences of SEDs shown in Figs. [5] 
and [7] (see also Tables [2] and |3]) . 



rameters of the jet in BL Lacs and quasars (Tables [5] and 
[3]). Most other parameter sequences produce correlations 
with opposite slopes than observed, while a variation of 
only the magnetic field over about two orders of magni- 
tude results in too small a range in radio-optical loudness 
to be consistent with the observed correlations (downward 
triangles in Fig. [5]). 

However, while reasonable agreement can be achieved 
with the R vs. L 7 correlation for quasars and BL Lacs and 
G r vs. L opt correlation for BL Lacs, our favored sequence 
of changing Lu (Lj) with eb = 1 predict steep-slope pos- 
itive correlations between G and Lvlba for quasars, in 
discrepancy with the slope of the observed positive correla- 
tion (Fig.jS] left panel). A better agreement is achieved for 
changes of the accretion power L d and the Lorentz factor 
of the jet (black circles and downward triangles in FigJTDJ). 
We therefore conclude that a simple one-parameter se- 
quence is not able to reproduce the observational results. 
The substantial scatter in all correlations also might be an 
indication that those are not single-parameter sequences. 
In a realistic scenario, we expect that a combination of 
different parameter variations and physical relations be- 
tween parameters are responsible for the observed range of 



radio, optical, and 7-ray luminosities. To disentangle the 
effect of multiple parameter sequences it is essential to pro- 
vide information on individual parameters (e.g. variance 
in Lorentz factor, variance in B field, etc.) from further 
studies. Moreover, the dispersion of the observed correla- 
tions can be used to constrain the distribution of initial 
parameters of the model. For this purpose, (quasi-) simul- 
taneous observations of the SEP (from radio to 7-ray) of 



more than 100 blazars (Planck Collaboration et al. 2011 



Giommi et al.ll201ll ) may be used to test the found corre- 
lations and to constrain further the jet model. This will 
be presented in a subsequent paper. 

Predictions of the model sequences are presented for 
quasars in Fig. (top panel) and BL Lacs in Fig. [TU] 
(empty symbols) . These can be tested with significant L opt 
- G r and Lvlba _ G correlations derived from larger 
samples of blazars. 



5. Summary 

Using the sample of 78 core dominated blazars (59 quasars 
and 19 BL Lacs) detected by Fermi /LAT from the 
MOJAVE-1 sample, we investigate the relations between 
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Fig. 9 Model correlations of the 7-ray-to-radio flux ratio vs. optical luminosity along sequences varying individual 
parameters. Each sequence represents the effect of a variation of one individual base parameter. The legend in the 
top left indicates the parameters varied along the sequences. Filled symbols indicate SEDs characteristic of FSRQs, 
while open symbols indicate sequences of SEDs characteristic of BL Lacs. The long-dashed line marks the power-law 
correlation that best fits the observational data. The sequences marked in blue dashed line (squares) correspond to 
the sequences of SEDs shown in Figs. [5] and [7] (see also Tables [H and [3]) . 



their radio, optical, and 7-ray emission available from non- 
simultaneous observations. Our main results are summa- 
rized as follows: 

— The detection rate of the 7-ray blazars from the radio- 
selected MO J AVE- 1 sample is high for radio bright 
blazars and it is independent of their optical nuclear 
fluxes. Most of non-detected 7-ray blazars are radio 
weak quasars. They are optically faint but also very 
luminous, L opt > IO 44 5 erg s _1 . 

— The known correlation between 7-ray flux and compact 
radio flux density (measured quasi-simultaneously) is 
also valid for non-simultaneous measurements of L 7 
and Lvlba for 59 quasars and 19 BL Lacs. This in- 
dicates that the Doppler-factors of the jet-plasma in 
the gamma and radio regimes are fairly constant on 
a timescale of a few years. The L 7 - Lvlba corre- 
lation is significantly stronger than that in the L 7 - 
Lopt plane. The latter holds exclusively for quasars 
at the confidence level (c.l.) of 98.9%. The relation 
logLvLBA = (0.83 ± 0.03) logL 7 + (5.96 ±0.26) can be 
used as a predictor of 7-ray flux from radio flux density 
at 15 GHz. 



— We report a statistically significant positive correla- 
tion (> 99%) between 7-ray luminosity and radio- 
optical loudness (R oc SvLBA/Sopt) f° r both quasars 
and BL Lacs. 

— Radio-optical loudness and 7-optical loudness of the 
population of quasars is one order of magnitude larger 
than those of BL Lacs, while no difference is found 
between 7-radio loudness of two populations. 

— We find a correlation between radio luminos- 
ity at 15 GHz and 7-ray-optical loudness (G oc 
S 7 /S opt ) for quasars, L V lba oc G°- 45±012 (c.l. > 
99.99%). The relationship between 7-ray-radio loud- 
ness (G r oc S^/SVlba) and optical nuclear lumi- 
nosity is marginally significant for BL Lacs (£ op t oc 
g-o.8i±o.38. c 1 gg g o^-j w jth a large error in the power 

index. 

— We employ the generic leptonic model for blazars 
to reproduce the observed correlations between 7-ray 
emission, optical, and radio emission described above. 
Among all parameters of the jet model, only a vari- 
ation of the accretion power (and the jet power) - 
with the magnetic field in equipartition with relativis- 
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Fig. 10 Model correlations of the 7-ray-to-optical flux ratio vs. radio luminosity along sequences varying Lp (for 
quasars, with fj — const.) and Lj (for BL Lacs), keeping the magnetic field at equipartition (squares), and varying B 
(downward triangles). Filled symbols indicate SEDs characteristic of FSRQs, while open symbols indicate sequences 
of SEDs characteristic of BL Lacs. The sequences marked in blue dashed line (squares) correspond to the sequences 
of SEDs shown in Figs. [B] and [7] 



tic electrons in the jet — can reasonably well repro- 
duce the observed strong correlations between L 7 and 
R for quasars and BL Lacs, L opt and G r for BL Lacs, 
and -Lvlba and G for quasars. We conclude that the 
change of the accretion (and jet) power is the main 
parameter responsible for the trends in most of the 
observed correlations but a combination of variations 
of the accretion power with other jet parameters such 
as the jet viewing angle, Lorentz factor, and the mag- 
netic field is needed to fit all correlations. Independent 
variations of the latter parameters will open up too 
large a parameter space, which makes infeasible the 
consideration of all possible combinations of parame- 
ter variations in this work. A larger sample of blazars, 
which will include the faint end of the population, is 
needed to better constrain the correlations between the 
bands. Ideally, one would need a simultaneous multi- 
band monitoring of a larger sample of blazars to test 
different jet models. 
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Table 4. Optical, radio, and 7-ray parameters of the 76 AGN from the MOJAVE-1 sample. Column (1) is the IAU 1950.0 
name of the object; (2) is the optical spectral type: quasars (Q), BL Lacs (BL), and radio galaxies (G); (3) is the redshift; (4) 
is the logarithm of the total radio flux (radio flux density multiplied by 15 GHz) in units of erg s _1 cm~ 2 ; (5) is the logarithm 
of the nuclear optical flux (optical flux density multiplied by 5100 A) in units of erg s _1 cm~ 2 ; (6) is the logarithm of the 7-ray 
energy flux between 0.1-100 GeV in units of erg s _1 cm~ 2 ; (7) is the logarithm of the radio luminosity at 15 GHz in units of 
erg s _1 ; (8) is the logarithm of the nuclear optical luminosity at 5100 A in units of erg s _1 ; (9) is the logarithm of the 7-ray 
luminosity in units of erg s _1 . All fluxes are Jf-corrected. 



Name 


Sp. type 


z 


log SyLBA 


!°g S'opt 


log S' 7 


log ivLBA 


log L opt 


log L 1 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


0059+581 


Q 


0.644 


-12.29 


-12.67 


-10.27 


44.96 


44.58 


46.97 


0106+013 


Q 


2.099 


-12.34 


-12.89 


-9.82 


46.19 


45.63 


48.70 


0109+224 


BL 


0.265 


-12.81 


-11.42 


-10.34 


43.52 


44.91 


45.99 


0133+476 


Q 


0.859 


-12.12 


-12.19 


-9.89 


45.44 


45.37 


47.67 


0202+149 


Q 


0.405 


-12.51 


-13.82 


-11.13 


44.25 


42.94 


45.63 


0202+319 


Q 


1.466 


-12.46 


-12.02 


-10.68 


45.68 


46.12 


47.46 


0212+735 


Q 


2.367 


-12.38 


-13.29 


-10.31 


46.27 


45.37 


48.34 


0215+015 


Q 


1.715 


-12.84 


-13.33 


-10.18 


45.46 


44.98 


48.13 


0234+285 


Q 


1.207 


-12.20 


-12.83 


-10.04 


45.72 


45.09 


47.89 


0235+164 


BL 


0.940 


-12.61 


-13.20 


-9.47 


45.05 


44.46 


48.19 


0316+413 


G 


0.018 


-11.80 


-10.35 


-9.76 


42.06 


43.51 


44.09 


0336-019 


Q 


0.852 


-12.46 


-12.14 


-10.70 


45.09 


45.41 


46.84 


0403-132 


Q 


0.571 


-12.67 


-12.24 


-10.88 


44.45 


44.88 


46.24 


0415+379 


G 


0.049 


-12.05 


-12.60 


-10.99 


42.69 


42.14 


43.76 


0420-014 


Q 


0.914 


-11.79 


-12.29 


-10.05 


45.84 


45.34 


47.57 


0422+004 


BL 


0.310 


-12.57 


-11.76 


-10.85 


43.92 


44.73 


45.63 


0458-020 


Q 


2.286 


-12.43 


-13.01 


-10.53 


46.18 


45.61 


48.08 


0528+134 


Q 


2.070 


-11.91 


-13.25 


-9.91 


46.60 


45.26 


48.60 


0529+075 


Q 


1.254 


-12.69 


-13.04 


-10.24 


45.28 


44.93 


47.73 


0529+483 


Q 


1.162 


-12.80 


-13.30 


-10.59 


45.09 


44.59 


47.29 


0605-085 


Q 


0.872 


-12.54 


-12.73 


-10.59 


45.04 


44.85 


46.98 


0736+017 


Q 


0.191 


-12.56 


-11.95 


-10.33 


43.45 


44.06 


45.68 


0748+126 


Q 


0.889 


-12.32 


-12.25 


-10.62 


45.28 


45.35 


46.97 


0754+100 


BL 


0.266 


-12.55 


-11.32 


-10.61 


43.78 


45.01 


45.72 


0805-077 


Q 


1.837 


-12.62 


-12.55 


-9.98 


45.76 


45.83 


48.40 


0808+019 


BL 


1.148 


-12.69 


-12.49 


-10.78 


45.19 


45.38 


47.09 


0814+425 


BL 


0.245 


-12.71 


-12.96 


-10.06 


43.55 


43.29 


46.19 


0823+033 


BL 


0.506 


-12.67 


-11.47 


-10.97 


44.32 


45.52 


46.02 


0827+243 


Q 


0.940 


-12.51 


-12.09 


-10.39 


45.14 


45.56 


47.27 


0829+046 


BL 


0.174 


-12.68 


-11.06 


-10.54 


43.24 


44.86 


45.38 


0836+710 


Q 


2.218 


-12.52 


-12.22 


-10.12 


46.07 


46.36 


48.46 


0838+133 


Q 


0.681 


-12.94 


-12.32 


-10.92 


44.36 


44.99 


46.39 


0851+202 


BL 


0.306 


-12.20 


-11.07 


-10.42 


44.28 


45.41 


46.06 


0906+015 


Q 


1.024 


-12.38 


-12.02 


-10.11 


45.37 


45.73 


47.64 


0917+624 


Q 


1.446 


-12.86 


-13.19 


-10.58 


45.26 


44.93 


47.54 


1055+018 


Q 


0.890 


-12.09 


-12.39 


-10.01 


45.51 


45.20 


47.59 


1124-186 


Q 


1.048 


-12.36 


-12.95 


-10.40 


45.41 


44.82 


47.38 


1127-145 


Q 


1.184 


-12.28 


-12.05 


-10.08 


45.62 


45.86 


47.82 


1156+295 


Q 


0.729 


-12.30 


-12.02 


-10.13 


45.09 


45.36 


47.26 


1219+044 


Q 


0.965 


-12.84 


-12.19 


-10.60 


44.84 


45.49 


47.09 


1222+216 


Q 


0.432 


-12.86 


-12.25 


-10.42 


43.96 


44.58 


46.41 


1226+023 


Q 


0.158 


-11.35 


-11.71 


-9.64 


44.48 


44.11 


46.19 


1228+126 


G 


0.004 


-12.39 


-9.78 


-10.76 


40.14 


42.76 


41.78 


1253-055 


Q 


0.536 


-11.55 


-11.65 


-9.36 


45.51 


45.40 


47.69 


1308+326 


Q 


0.997 


-12.34 


-13.09 


-10.03 


45.38 


44.63 


47.69 


1324+224 


Q 


1.400 


-13.02 


-12.51 


-10.43 


45.07 


45.58 


47.66 


1334-127 


Q 


0.539 


-11.86 


-12.59 


-10.38 


45.19 


44.47 


46.68 
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Table 4. continued. 



Name 


Sp. type 


z 








log LvLBA 


log L opt 


logL 7 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1502+106 


Q 


1.839 


-12.75 


-13.08 


-9.08 


45.63 


45.30 


49.30 


1510-089 


Q 


0.360 


-12.35 


-11.61 


-9.15 


44.29 


45.03 


47.49 


1546+027 


Q 


0.414 


-12.36 


-12.27 


-10.46 


44.42 


44.51 


46.33 


1548+056 


Q 


1.422 


-12.35 


-12.58 


-11.01 


45.75 


45.53 


47.10 


1606+106 


Q 


1.226 


-12.46 


-12.52 


-10.41 


45.48 


45.42 


47.54 


1611+343 


Q 


1.397 


-12.16 


-12.41 


-11.09 


45.93 


45.67 


46.99 


1633+382 


Q 


1.814 


-12.18 


-12.66 


-9.86 


46.19 


45.71 


48.51 


1641+399 


Q 


0.593 


-11.87 


-11.66 


-10.02 


45.29 


45.50 


47.14 


1726+455 


Q 


0.717 


-12.47 


-12.66 


-10.58 


44.89 


44.70 


46.78 


1730-130 


Q 


0.902 


-11.77 


-13.55 


-10.37 


45.84 


44.06 


47.24 


1749+096 


BL 


0.322 


-12.07 


-12.41 


-10.10 


44.46 


44.12 


46.42 


1803+784 


BL 


0.680 


-12.42 


-11.45 


-10.36 


44.88 


45.86 


46.95 


1807+698 


BL 


0.051 


-12.68 


-11.84 


-10.50 


42.10 


42.94 


44.28 


1823+568 


BL 


0.664 


-12.45 


-12.58 


-10.41 


44.83 


44.71 


46.87 


1828+487 


Q 


0.692 


-12.70 


-12.08 


-11.12 


44.63 


45.24 


46.21 


1849+670 


Q 


0.657 


-12.57 


-11.73 


-9.78 


44.70 


45.54 


47.49 


1958-179 


Q 


0.650 


-12.39 


-12.29 


-10.29 


44.87 


44.96 


46.96 


2131-021 


BL 


1.285 


-12.47 


-12.76 


-11.08 


45.52 


45.24 


46.91 


2145+067 


Q 


0.990 


-11.80 


-11.91 


-10.81 


45.91 


45.80 


46.90 


2155-152 


Q 


0.672 


-12.48 


-12.25 


-10.89 


44.81 


45.04 


46.41 


2200+420 


BL 


0.069 


-12.06 


-11.29 


-10.09 


42.99 


43.76 


44.96 


2201+171 


Q 


1.076 


-12.51 


-12.98 


-10.16 


45.29 


44.82 


47.64 


2209+236 


Q 


1.125 


-12.61 


-13.09 


-10.88 


45.24 


44.76 


46.97 


2223-052 


Q 


1.404 


-12.13 


-12.52 


-10.31 


45.96 


45.57 


47.78 


2227-088 


Q 


1.560 


-12.48 


-12.46 


-9.83 


45.72 


45.75 


48.38 


2230+114 


Q 


1.037 


-12.15 


-12.13 


-10.04 


45.61 


45.64 


47.72 


2251+158 


Q 


0.859 


-11.81 


-11.88 


-9.09 


45.75 


45.68 


48.47 


2331+073 


Q 


0.401 


-12.72 


-11.73 


-10.82 


44.03 


45.02 


45.93 


2345-167 


Q 


0.576 


-12.41 


-12.53 


-10.63 


44.72 


44.60 


46.50 



